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Abstract 

In a recent paper, a minimal supersymmetric (SUSY) 50(10) x ^4 based unified model of 
flavor for quarks and leptons was proposed with two 10 and one 126 contributing to fermion 
masses. An important aspect of this model is that Yukawa couplings emerge dynamically from 
minimization of the flavon potential, thereby reducing the number of parameters considerably. 
We make a detailed numerical analysis of this model for fermion mixings including SUSY thresh- 
old effects at the TeV scale and type-I corrections to a type-II dominant seesaw for neutrino 
masses. This is a single-step breaking model with SUSY SO{10) broken at the Grand Unified 
Theory (GUT)-scale of 2 x 10^*^ GeV to the Minimal Supersymmetric Standard Model (MSSM). 
The model has only 11 parameters, and therefore, the charged fermion fits predict the masses 
(up to an overall scale) and mixings in the neutrino sector. We present correlations for the 
different predictions in the neutrino mixing parameters. Typically, we find ^13 ~ 5°, Ou ~ 32° 
and 023 > 48°. We also find that proton decay mode p — )■ K'^u^ has a partial lifetime of ~ 10^^ 
yrs, which is within reach of the next round of planned proton decay searches. The successful 
fit for fermion masses requires the Higgs mass to be below 129 GeV in this model. If the Higgs 
mass lies between 120-128 GeV, as suggested by the recent LHC data, we find a lower limit on 
the light stop mass of 755 (211) GeV for /x > (< 0). 
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I. INTRODUCTION 



The measurement of neutrino masses and mixings during the past decade has 
provided the first evidence for physics beyond the standard model (SM). It has also 
raised the possibility that this new knowledge may help to unravel the physics of fiavor 
and also possibly unlock the mystery behind the origin of matter in the universe. An 
important recent experimental finding in this area has been the announcements by the 
T2K [li, MINOS Q and Double CHOOZ [s] experiments that one of the hitherto unknown 
neutrino mixing angles, namely 6'i3, may not only be non-zero but could be "large". The 
result is, however, still preliminary, and more results are awaited from ongoing as well as 



near future experiments jj] . Some hints for a large 6^13 have also been suggested by a global 
analysis of the existing oscillation data Is, A non-zero 613 has profound implications 
for our understanding of the physics of neutrino mass, and it is therefore timely to search 
for the prediction of ^13 in various models 3]- 

A key question before theorists now is: what is the big picture of fiavor where 
known quark and lepton masses and mixings fit together? A framework that suggests 
itself is grand unified theories (GUT), where all matter and all forces (except gravity) 
unify at high energy scale. Since single-step coupling unification in these theories requires 
supersymmetry (SUSY), we will focus on SUSY-GUTs for approaching the fiavor problem 
and use the seesaw mechanism Sj to understand the small neutrino masses. A possible 
advantage of this framework is that matter unification is expected to reduce the number 
of free parameters that determine fermion masses and mixings from 31 parameters in 
the seesaw extended standard model, so that one can make predictions to make the 
model testable. The minimal scenario which is predictive for neutrinos is a renormalizable 
supersymmetric 50(10) model with 10 and 126 Higgs fields_^ contributing to fermion 
masses [9]. This model embodies both the type-I and type-! 
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lOj seesaw contributions to 
Fitting charged fermion 



neutrino masses and has been analyzed in many papers 
masses in these models leads to predictions for lepton mixing angles and will be tested with 
higher precision measurement of these angles and the neutrino mass hierarchy. Detailed 
analysis of the superpotential and symmetry breaking for these models have been carried 
out in Ref. Q. 
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The next step with these models is to see if Yukawa couphngs can be predicted 
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as consequences of higher-scale symmetries. There have been several such approaches 
[isl l which adopt the point of view that Yukawa couplings are dynamical fields, namely 
the flavon fields, whose vacuum expectation values (vevs) determine the observed Yukawa 
couplings. If these vevs could be the result of the minimization of simple superpotentials, 
that would reduce the number of parameters and would indeed be a step further in the 
search for an understanding of flavor. In a recent paper |16i], such a model was presented 
where the GUT scale 5*0 (10) theory was extended to include three fiavon fields with an 
5*0(10) X 5*4 symmetry, where the fiavon fields, as well as the three matter families, form 
a three dimensional irreducible representation of 84^. 

It was shown in [I6!] that the ground state of the 5*0(10) x 54 theory has only 11 
parameters describing all flavor, i.e. quark and charged lepton masses and quark mixings 
as well as neutrino masses and lepton mixings. Using pure type II seesaw contribution, 
it was argued that the model is in qualitative agreement with observations. In this paper 
we present a detailed numerical analysis of this model and its predictions in the neutrino 
sector as well as for proton decay. We find that once we include the SUSY threshold 
corrections to quark masses and a small type-I contribution to neutrino masses, all existing 
data in the quark and charged lepton sector can be fitted to a good accuracy; furthermore, 
the model predicts the solar to atmospheric mass ratio AitLq/ Am'^^^ and the solar and 
atmospheric mixing angles in agreement with the current neutrino oscillation data. In 
addition, it predicts the reactor mixing angle 6*13 as well as the Dirac and Majorana phases 
in the leptonic mixing matrix, which can be used to test the model. In particular, some 
outstanding features of the predictions are: a value of 623 > 48°, 4.5° < ^^13 < 5.4°, and 
a Dirac CP phase of 144°. We will demonstrate the determination of these features, and 
also note the correlation between the different parameters, which can make it easier to 
rule out the model. 

We also find that in order to get the desired threshold corrections to the b- 
quark mass to fit observations, we need a large negative fi and/or A-terms in the model. 
We discuss predictions for Higgs mass for this choice of MSSM parameters, and find 
that the Higgs mass should be below 129 GeV in order to satisfy the constraints from 
the fermion sector fit. We also obtain lower limits on the squark masses from the same 



constraints. In particular, we note that if the Higgs mass is discovered between 120-128 
GeV, the hght stop should be heavier than 755 (211) GeV for /i > (/x < 0) in this model. 
These features could be used to test the model at the LHC. 

Finally, we discuss proton lifetime predictions in this model. We find that in 
order to suppress the RRRR contributions to proton lifetime, we need to work in the low 
tan /3 regime, which is the reason we need a large A-term for the threshold correction, 
as noted above. We find that the S-violating LLLL - terms contribute dominantly to 
p — > K^v^ decay mode for which we find a partial lifetime of ~ 10'^^ yrs. We emphasize 
that proton lifetime predictions as well as the predictions for neutrino mixings can be 
used to test the model. 

The paper is organized as follows: in Section II, we present the essential points 
of the minimal 5*0 (10) x model; in Section III, we discuss the fermion mass fits and 
the predictions for the neutrino sector; in Section IV, we discuss the SUSY threshold 
correction required to fit quark masses and its implications for gluino, stop and the Higgs 
masses in the model. In Section V, we discuss our predictions for proton lifetime and in 
Section VI, we summarize our results. 

II. DETAILS OF THE MODEL 

The class of 5*0(10) models we are interested in here have two 10 Higgs fields 
(denoted by H,H') and one 126 (accompanied by 126, denoted by A + A). The 5*0(10)- 
invariant Yukawa couplings of the model are given by: 

Cy = hijjipH + h'i)i)H' + f^^A (1) 

where ^/''s denote the 16 dimensional spinors of 5*0(10), which contain all the matter fields 
of each generation; there are, of course, three such fields, though we have suppressed the 
generation indices. The Yukawa couplings are therefore 3x3 matrices in generation space. 

The effective Yukawa couplings /, /i, h' are assumed to have descended from a 

higher scale theory which has an 5*4 symmetry broken by flavon fields with particularly 

[ I 

aligned vevs (see e.g. Ref. [16?]). We do not need to know the detailed form of these flavon 
interactions for our analysis in this paper, and we will simply write down the effective 
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form of h,h' , and / that follow from it. Before doing that, we wish to point out that it 
is the / coupling which is responsible for neutrino masses via seesaw mechanism; it also 
contributes to charged fermion masses. 

The fermion mass matrices are derived from the Yukawa interaction as fol- 
lows: after GUT symmetry breaking, two linear combinations of the SM doublets remain 
light, denoted by and Hd- Typically, = J2aUu{d)aHu(d)a, where Hu(d)a are the 

up(down)-type SM doublets in the {H,H',A). The effective Hu coupling at the MSSM 
scale is then given by QHuu'^{hUuH + /C4a + h'UuH'), and similarly for Ha- The fermion 
mass matrices can be written in terms of these couplings as 



Mu 
Md 

Me 



tan P 
ri 



tan (3 



(2) 



where we have absorbed the mixings Uda) and the vevs {Hu,d) = K'u^-i with tan/3 
i^u/i^d-i to re-define the Yukawa coupling matrices as follows: 



with the ratios 



h = KuUuHh, / = /, h' = h , 

n ri 



UdH UuA UuH' 

UuH t^dA 'JdH' 



(3) 



(4) 



These effective coupling matrices determined by the flavon sector are of the 



form 



h 






M 



/ 



(5) 



mi mi 
mi mo mi — mo 

y mi mi — mo mo j 



(6) 
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^0 5-6^ 

h' = 6 . (7) 
y-5 J 

It was shown in Ref. [16| that 5*4 symmetry constrains h to have the above rank-one form. 
The parameters mQ,mi,6 are chosen to be complex, giving a total of 10 parameters in 
the charged-fermion sector. 

The neutrino mass matrix is, in general, given by a combination of the type-I 
and II seesaw mechanism: 



(8) 



where vl^r are the vevs of the SM triplet Higgs fields in 126. If we assume type-II 
dominance, imposed by the ratio of vl and vr, and the magnitude of the coupling /, the 
neutrino mass matrix takes an approximate tri-bimaximal (TBM) Q form via the 
form of / given in Eq. ([6]), which is diagonalized by 

/ /I /I ^ 



TBM 



V 



(9) 



Note that the full neutrino mixing matrix, U pmns = VeVu, where Vi and are the unitary 
matrices that diagonalize the charged lepton and neutrino mass matrices, respectively. 
Hence, we will necessarily have corrections to the TBM mixing given by Eq. coming 
from the charged lepton mass matrix as well as the type-I contribution. Note further 
that since the / matrix also contributes to the quark and charged lepton masses, neutrino 
masses and quark masses are connected, thus making the model predictive. 

III. PREDICTIONS IN THE NEUTRINO SECTOR 



We diagonalize the mass matrices given by Eqs. ([2]) and use the Minuit2 tool 



library 20| to minimize the sum of chi-squares for the mass eigenvalues and CKM mixing 
in the charged-fermion sector as well as the mass-squared differences AitLq and Aml^^ 
in the neutrino sector. We note that a small but non-zero type-I contribution is required 



in the neutrino mass matrix given by Eq. (|8]), in order to have a consistent fit with the 
correct mass squared ratio AmQ/Am^^.jj^; on the other hand, a large type-I contribution 
will spoil the TBM structure given by Eq. ([6]) and hence results in too small mixing angles. 
A balancing of the two is required in order to satisfy the observed neutrino oscillation 
data. We also note that SUSY threshold corrections to the down quark mass matrix 22 1 
must be included in order to get a consistent fit for the charged fermion sector. The 
details of this analysis are given in Section IV. Also, as discussed in Section V, the proton 
decay constraints are satisfied only for low tan /3 in this model, as the Yukawa couplings 
responsible for the proton decay rates grow with tan /3. Hence, we have chosen tan /3 = 10 
for our numerical analysis given below. 



M (GeV) 


83.06 


mo (GeV) 


1.201 - 0.9007i 


mi (GeV) 


0.2033 - 0.01170i 


5 (GeV) 


0.2129 + 0.0820n 


ri/ tan f3 


0.01624 


r2 


-0.1382 


rs 


0.1358 


a 


5.0° 



TABLE I. Best fit values for the model parameters at the GUT scale. 



The fit results are displayed in Tables I and II; Table I gives the numerical 
values of the model parameters yielding the best fit values shown in Table II. Here, a is 
the mixing angle for the third generation matter fermion ip with the vector-like field ipv 
specific to the model [l6|; the limit a = gives the form for the mass matrices dictated 
by 5*4 symmetry, as given by Eqs. (E]) exactly, and the fit value of a = 5° approximates 
this limit. With this in mind, note that the top quark mass in the model is given by 

mt ~ UuHK'uhss cos^ a (10) 

In the neutrino sector, as noted earlier, the correct mass squared ratio 
Am^/Am^^j^ as well as large solar and atmospheric mixing angles fix the relative size 



me (MeV) 
(MeV) 
(GeV) 
TTid (MeV) 
rus (MeV) 
rub (GeV) 
m^, (MeV) 
mc (MeV) 
mt (GeV) 

Kb 

Vcb 

J 

Aml/Aml^^ 



best fit 



0.3585 
75.6717 
1.2922 
2.0034 
23.4494 
1.0335 
0.8192 
207.4990 
82.8964 
0.2245 
0.0034 
0.0351 
2.052 X 10-5 
0.0311 



3.39 



exp value 



0. 35851°:r 

1 9Q99+0.0013 

1. zyzz_Q QQ^2 

1.5036l°;ig^ 

i.UD,3D_Q 0865 

0.7238l[J:li^ 
210.3273;i?:0036 

Qo /| 000+30.2676 
"^•^"J-^ "5-14.7686 

0.2243 ± 0.0016 
0.0032 ± 0.0005 
0.0351 ±0.0013 
(2.2 ±0.6) X 10-5 
0.0320 ± 0.0025 



TABLE II. Best fit values for the charged fermion masses and the most relevant quark mixing 
parameters, as well as the solar-to-atmospheric mass squared ratio. The la experimental val- 
ues 



231 ] ■ with masses and mixings extrapolated by MSSM renormalization group (RG) equations 
to the GUT scale, are also shown for comparison. Note that the values of the bottom quark 
mass and the CKM mixing parameters involving the third generation quoted here include the 



SUSY-threshold corrections (see Section IV). 



between the type-I and type-II contributions, and then the overall scale is determined 
from the largest mass eigenvalue, assuming a normal hierarchy. We find that for the best 
fit parameters shown in Table I, and for vr = 2.0 x 10^^ GeV (same as the GUT scale), 
vl = 6.810 eV yields the right neutrino mass scale with ~ 0.05 eV. For estimating the 
proton decay rates as well as for the neutrino masses and mixing, we must extract the 
magnitudes of the raw yukawa couplings h, f, h', which can be done using the expressions 
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in Eq. (|3]). However, these couplings depend on the vev mixing parameters Uqa, and 
hence, there is some freedom in their determination, although the unitarity constraints 



on the f/'s, J2a l^<?«P ^ 1; and the top-quark mass relation in this model, given by 
Eq. ( JTOl) . provide some restriction on these mixing parameters. The values chosen for the 
up-type mixings are Uuh = 0.40, Uua = 0.4033, and UuH' = 0.72, and using the fit values 
for the r's from Table 1 and Eq. (jlj), the resulting values for the down-type mixings are 
UdH = 0.06497, f/(iA = —0.4739, and UdH' = 0.8611. Given these values and the running 
vevs = 123.8, Kd = 17.9 GeV at GUT-scale [231], the resulting dimensionless couplings 
are found to be 



h 











V 



1.677 



/ 



-5.628 + 0.3238i) x 10"^ 
;-3.326 + 2.494z) x lO^^ 
(2.763 - 2.461i) x lO^^ 



-5.628 + 0.3238i) x 10"^ ^ 
(2.763 - 2.461i) x lO^^ 
;-3.326 + 2.494i) x lO^^ J 



h' 



I 



( 

(-5.628 + 0.3238?) x 10"^ 
\^ (-5.628 + 0.3238i) x 10"^ 

( (3.243 + 1.250i) x 10"^ (-3.243 - 1.250i) x 10"^ ^ 

(3.243 + 1.250i) x 10"^ 

\^ (-3.243 - 1.250i) X 10"^ 

The predicted values for the neutrino mixing parameters corresponding to the 
best fit parameter values in the model given in Table I are summarized in Table III; 
the correlations between the different parameters in the neutrino sector are shown in the 
scatter plots of Figure [H We find that a consistent fermion sector fit in this model predicts 
the reactor mixing angle ^13 to be non-zero and within a very narrow range 4.5° — 5.5°, 
which, though not "large" , is well within the 3(T lower bound of the recent experimental 
results [l|-|3|. The model also predicts a large atmospheric mixing angle, ^23? and a specific 
correlation with the solar mixing angle (Fig. [1]). These predictions can be used to test 
the model in the ongoing and near future neutrino experiments. 



Ill 
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predicted value 


3a exp range 


^12 


32.34° 


(30.6 - 36.8)° 




49.41° 


(35.7 - 53.1)° 


^13 


5.13° 


(1.8 - 12.1)° 


Sd 


144.4° 





TABLE III. The model predictions for the neutrino mixing angles for the best fit parameter 
values given in Table I. We also show the 3a range of values from the updated global neutrino 
data analysis [sj. 
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FIG. 1. The predicted correlation between the neutrino sector observables, namely Attiq/ Arn^^^ 
and the mixing angles ^12)^235^135 based on the fermion sector fit in the model. The 1,2 and 
3a experimental limits for ^12, ^23 s>.nd AuIq/ Arn^^^^ are also shown (shaded regions). 
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IV. THRESHOLD CORRECTION AND LOW ENERGY MASS SPECTRUM 



In order to compare the fermion masses and mixing values obtained from the 
model at the GUT scale with the experimental values at the weak scale, we must take into 



21^. 



22 1 . There are two main contributions 



account the SUSY threshold correction effects 
to the SUSY threshold correction to the fermion masses: one coming from the gluino loop 
and another from the chargino loop. The largest correction is for the bottom mass, which 



is given by 



2l| 



1 

16^ 

and the function is given by 



^^ei + e2 Vtb ', 12 
where ei = — ^/xm^ tan/3 J3(m?, mf^, m|), (13) 



62 = --^^iAtyl tan /3 m| , m|) (14) 



aMog (f) +&C log (^)+ca log (^) 
^' = {a-h)ih-c){a-0 (^'^ 

Similarly, if we do not add any off-diagonal threshold corrections, the CKM mixings 



involving the third generation receive corrections as follows 21| : 

5Vub SVcb SVtd SVts 



62 (16) 



Vub V,b Vu Vts 
However, once off diagonal threshold corrections are included, there are further changes 
to CKM mixing, which we take into account in our numerical analysis. 

From the numerical fit, we find that at the GUT scale, without including the 
threshold corrections, some of the best-fit values predicted by the model do not agree with 
experimental values extrapolated to the GUT scale (see Table II). In particular, we find 

rUb = 1.37 GeV, |Kfe| = 0.0015, = 0.0160, \Vtd\ = 0.0047, \Vts\ = 0.0153. (17) 

Comparing these values with the experimental values, we note that large negative thresh- 
old corrections are required for the model to have a consistent fermion-sector fit. We 
parameterize the SUSY threshold corrections at the GUT scale by modifying the third 



12 



generation elements of the down-quark mass matrix as follows: 

^0 5i3^ 

623 (18) 

V 5l3 S23 §33 J 

and ^Ad is given by Eq. ([2]). The required threshold corrections at the GUT scale are: 



M'a = Md + 5Md, where SMc 



5i3 = 0.09 GeV, 623 = -0.96 GeV, ^33 = -20.68 GeV. 



(19) 



Note that these threshold corrections, when extrapolated down to the weak scale and 
added to the RG-extrapolated values of the 6-quark mass and the CKM mixing parame- 
ters, yield results within la range of the experimental values at Mz- 

At the weak scale, it is clear from Eqs. (fT3|) and (fl^ that for the large negative 
threshold corrections given by Eq. flTIJl) . we must have fi < if the gluino term is dominant, 
or opposite signs for /x and At if the chargino contribution is dominant. These observations 
have important consequences for the MSSM light Higgs mass as well as the sparticle 
spectrum, as shown below. 

The one-loop radiative correction to the MSSM light Higgs mass is given by 



25| 



ml ~ 



1 + 



rriz 



1 + 



(20) 



where m^^ ^ the stop mass eigenvalues, which are obtained by diagonalizing the matrix 
/ 



Ml 



mf^ +mf + m\ cos 2/3 (| - | sin^ 6 
At + fi cot (3) 



w 



mt[ 



■mt{—At + jjL cot (3) 
mf^ + mj + m| cos 2/3 (^| sin^ 6 



w , 



Note that the same mass parameters, namely the gluino and third generation 
squark masses as well as the trilinear term At, appear in the threshold correction and 
Higgs mass correction. Given that this model requires a large threshold correction to 
have a consistent fit in the fermion sector, we expect some correlation between the two 
corrections. To show it quantitatively, we have chosen the simplest case of mSUGRA- 
type GUT scale spectrum, as an illustration, although our results do not depend on the 
assumption of mSUGRA. Here we scan over the parameter space for 200 GeV < niQ < 
2 TeV, 200 GeV < mi/2 < 2 TeV, -5 TeV < Aq < 5 TeV and tan/3 = 10, using the 
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ISA JET package 26|]. The results are shown in Fig. [21 from which we find that the hghtest 
Higgs mass is required to be below 129 (128) GeV for < (> 0) to satisfy the bottom- 
quark mass constraint (vertical red shaded region) in this model. The horizontal green 



shaded region shows the range o 
recently reported at the LHC 



27]. 



Higgs mass in which a mild excess of events has been 



130 
128 



124 
122 
120 
118 
116 




FIG. 2. Higgs mass prediction in the model required by the bottom quark mass fit, including the 
threshold correction effects. The dashed vertical line is the best fit value for nrn,{Mz) without 
the threshold corrections. The vertical (red) shaded region is the la experimental range of 



mb{Mz) = (2.85 ±0.32) GeV, obtained from the LEP data [2j]. The horizontal (green) shaded 
region shows the range in which ~ 2>a excess of events for SM-like Higgs were observed recently 
at the LHC 

The correlations between the gluino mass and the light stop and sbottom 
masses for the required threshold corrections (shaded region in Fig. [21) are shown in 
Fig. [3l We find that the large threshold correction requirement in this model requires 
the gluino to be always heavier than the light stop, but not necessarily heavier than the 
light sbottom. Moreover, for gluino masses satisfying the current LHC lower bound of 



1.1 TeV |28|] and for Higgs mass between 120-128 GeV, we find a lower limit for the stop 
mass of 755 (211) GeV for yU > 0(< 0) and similarly for the sbottom mass of 1013 (895) 
GeV. The milder limit on the squark masses for /i < is because of the fact that in 
this case, the required negative threshold corrections can be obtained from both gluino 
and chargino contributions (c.f. Eq. ( lT3l) and ( lT4l) ). thus allowing for the At,Ah values 
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necessary for light stop and sbottom masses, respectively. However, for /i > case, the 
gluino contribution is of the wrong sign, and hence we must have very large negative At 
values to obtain the required threshold corrections. 



tl< 
U > 



3.5 - 

3 - 



.''Oh 





m; (TeV) 



me, (TeV) 



FIG. 3. The correlation between the gluino mass and the light stop and sbottom masses, 
corresponding to the shaded region in Fig. [2j The dashed line is for gluino mass equal to the 
relevant squark mass in the plot. The red shaded region is ruled out by the current LHC 
data [28 1 . 



V. EXPECTATIONS FOR PROTON DECAY 

We now turn to a discussion of proton lifetime in our model. As experimental 
limits on proton lifetimes keep increasing, many simple SUSY-GUT models have either 
been ruled out or become more and more constrained. It is therefore important to ensure 
that any GUT model for neutrino masses is consistent with those limits. As is well known, 
the dominant contributions to proton decay comes from color triplet Higgsino exchange 
in these modes [29|] and can in general lead to large proton decay amplitudes [30|. It 



was suggested in Ref. [3l|] that one way to suppress this amplitude without invoking 
cancellations is to choose appropriate flavor structure for the Yukawa couplings. The 
current model falls into this category where the existence of small elements in h, h' and 
/ matrices leads to the expectation that we should be able to satisfy the proton decay 
constraints without any cancellation. Let us now see how this occurs using the fit for 
Yukawa couplings we obtained in the previous section. 

15 



The colored Higgs triplets, (l)T + (t>T- ((3, 1, —1/3) + c.c.) responsible for proton 



decay in our model arise from 10 + 10' + 126 + 126 + 210 multiplets. Once the triplet 
fields and 0y are integrated out, both LLLL (C^) and RRRR {Cr) operators leading 
to proton decay emerge : 

W, = \cf'Euqim, + Cf'elu'iu'id^^ (21) 
The color triplet f fields are linear combinations of six fields, two of them arising from 



two lO's, three of them arising from 126 + 126 and one from 210. This leads to a 6 x 6 
dimensional mass matrix for the triplets: {4>f)a{MT)ab{4'T)b- One can write the dimension 
five operators in terms of the couplings h, h' and / as follows: 



C'i = chijhki + xihijh'i^i + X2h[jhki + x^K^h'^^i + Xifijfki 

+ x^fijhki + x^hijfki + Xfh'-jfki + Xsfijh'j^i- (22) 

Similarly, we can write Crrrr operator (and change Xj's to different coefficients Ui). The 
coefficient c is (M^ ^)ii and the coefficients Xj and Ui are also given by the components of 
M^^ . The proton decay amplitude can be written as 

^ ~ 4^74 ' ^^^^ 

where A = cAhh + xiAhj + ... using the c's and XjS given in Eq. fl22l) . and f3p is the nucleon 
matrix element of the three-quark operator. In our calculation we use rrig = 1.3 TeV as a 
typical value for the first two generation squark masses (to satisfy the LHC lower bound), 
and similarly m^r = 200 GeV as a typical wino mass. The proton decay width can be 
written simply as : 

r . ,2.T X 10- GeV,|Cn/(F.fl,P fiiiiT^V^ f fi^V .(24) 



Mt J V200 GeVy V ^9 
where f{F, D) are the hadronic form factors, typically of 0{1). Assuming that the colored 
Higgs from lO's are the lightest, we find the largest contribution to the p — )■ v^K^ to be 
arising from Using the values of h' from the Eq. ffTIl) and varying Mt within a factor 
5 of the GUT scale, we find that the partial lifetime of proton in this mode can be as large 
as 7 X 10^^ years. If however, we lower tan /9 to 5 {h'12 will be smaller in this case), or raise 
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the squark mass to 1.5 TeV, the hfetime can be be as large as 10^^ years. The partial 
lifetimes for other flavors of v are found to be larger. Similarly, the partial lifetimes for 
other decay modes, e.g., n — )■ z/^7r° are found to be much larger of (9(10^^) years. We 
should emphasize here that, in order to generate these numbers, we did not invoke any 
cancellation. The smallness of the elements of the Yukawa coupling matrices are sufficient 
in suppressing the decay rates. We note that the lifetime of the mode p — > v^K^ is within 
the search limit of the proposed Hyper-Kamiokande experiment 32|. 



VI. SUMMARY 



We have analyzed the predictions of a minimal 5*0(10) x 5*4 model of flavor 
with dominant type II seesaw form for neutrino masses, where the forms of the Yukawa 
couplings for the two 10 Higgs flelds and one 126 Higgs fleld are determined dynamically 
by flavon vevs at the minimum of the S'4-invariant flavon potential. The model has 
eleven parameters including complex phases and is, thus, a relatively economical one 
when compared to other models discussed in the literature. It gives a very good flt to the 
charged fermion masses and the CKM parameters, and it also predicts the neutrino mixing 
angles ^12, ^23 as well as Amg/Am^^j^ in agreement with observation. Furthermore, it 
predicts a non-zero value for ^13 between 4° — 6°, which is in the current experimentally 
preferred range. With more accurate determination of ^^13 - especially its correlation with 
^23, which our model predicts to be strictly larger than 45° - the model could be tested in 
near future. The model also predicts a normal hierarchy for the neutrinos and hence an 
effective neutrino mass in neutrino-less double beta decay which is a few milli-electron- 
volts and is thus not observable in the current round of the searches for this process. The 
proton lifetime for p — )■ Vy}z^ decay mode can be 10^"^ years. Finally, the successful flt for 
fermion masses require the Higgs mass to be below 129 GeV in this model and put lower 
bounds on the third generation squark masses which are well within the reach of LHC 



17 



VII. ACKNOWLEDGMENT 



The work of P. S. B. D., R. N. M., and M. S. is supported by the National 
Science Foundation grant number PHY-0968854 and the work of B. D. is supported by 
the DOE grant DE-FG02-95ER40917. P. S. B. D. and M. S. would like to thank Shabbar 
Raza and Mike Richman for their patient help with the numerical tools. 



[1] K. Abe et al [T2K Collaboration], Phys. Rev. Lett. 107, 041801 (2011) |arXiv: 1106. 2822] 
[hep-ex]]. 

[2] P. Adamson et al. [MINOS Collaboration], Phys. Rev. Lett. 107, 181802 (2011) 

[arXiv:1108.001 5' [hep-ex]]. 
[3] H. De. Kerrect [Double CHOOZ Collaboration], talk at the LowNu conference in Seoul, 

Korea (2011), |http : //workshop .kias .re .kr/lowiiull/?Program| 
[4] H. Steiner [Daya Bay Collaboration], Prog. Part. Nucl. Phys. 64, 342-345 (2010); J. K. 

Ahn et al. [RENO Collaboration], arXiv: 1003.139T] [hep-ex]: P. Novella (Double CHOOZ 

Collaboration], arXiv: 1105.6079 [hep-ex]. 
[5] G. L. Fogh et al, Phys. Rev. D 84, 053007 [arXiv:1106.6028] [hep-ph]]. 
[6] M. C. Gonzalez-Garcia, M. Maltoni, J. Salvado, JHEP 1004, 056 (2010) |arXiv: 1001. 4524) 

[hep-ph]]; T. Schwetz, M. Tortola and J. W. F. Valle, New J. Phys. 13, 109401 (2011) 

[arXiv:1108.1376i [hep-ph]]. 
[7] H. -J. He, F. -R. Yin, Phys. Rev. D 84, 033009 (2011) |arXiv:1104.2654 ' [hep-ph]]; Z.- 

z. Xing, Chin. Phys. C 36, 101 (2012) [ arXiv: 1106. 32441 [hep-ph]]; N. Qin and B. Q. 

Ma, Phys. Lett. B702, 143 (2011) [arXiv: 1106. 3284) [hep-ph]]; Y. -j. Zheng and B. Q. 

Ma, Eur. Phys. J. Plus 127, 7 (2012) [arXiv: 1106.4040) [hep-ph]]; E. Ma and D. Weg- 

man, Phys. Rev. Lett. 107, 061803 (2011) [arX iv:1106.4269l [hep-ph]]; X. -G. He and 

A. Zee, Phys. Rev. D 84, 053004 (2011) [arXiv:1106.4359l [hep-ph]]; S. Zhou, Phys. 

Lett. B704, 291 (2011) [a rXiv: 1106.4808) [hep-ph]]; T. Araki, Phys. Rev. D 84, 037301 

(2011) t^iv:1106.5211 [hep-ph]]; N. Haba, R. Takahashi, Phys. Lett. B702, 388 (2011) 

[arXiv:1106.5926) [hep-ph]]: D. Meloni, JHEP 1110, 010 (2011) [arXiv:1107.0 221 [hep-ph]]; 



18 



W. Chao, Y.-j. Zheng. [a?Xrv:1107.0738l [heo-ph]: H. Zhang, S. Zhou, Phys. Lett. B704, 
296 (2011) (arXiv: 1107. 10971 [hep-ph]]; X. Chu, M. Dhen and T. Hambye, JHEP 1111, 
106 (2011) [arXiv: 1107. 1589] [hep-ph]]; P. S. Bhupal Dev, R. N. Mohapatra, M. Sever- 
son, Phys. Rev. D 84, 053005 (2011) !arXiv:1107.2378J [hep-ph]]; R. d. A. Toorop, F. 
Ferugho, C. Hagedorn, Phys. Lett. B703, 447 (2011) [arXiv: 1107.3486 ^ [hep-ph]]; S. An- 
tusch, V. Maurer, Phys. Rev. D 84, 117301 (2011) [arXiv: 1107.3728) [hep-ph]]: W. Rode- 
johann, H. Zhang and S. Zhou, Nucl. Phys. B 855, 592 (2012) [ arXiv: 1107.39701 [hep-ph]]; 
Y. H. Ahn, H. -Y. Cheng and S. Oh, Phys. Rev. D 84, 113007 (2011) [ arXiv:1107.4549l 
[hep-ph]]; S. F. King and C. Luhn, JHEP 1109, 042 (2011) [arXiv: nOL5332| [hep-ph]]; Q. 
-H. Cao, S. Khahl, E. Ma, H. Okada, Phys. Rev. D 84, 071302 (2011) [a rXiv: 1108. 0570) 
[hep-ph]]; D. Marzocca, S. T. Petcov, A. Romanino and M. Spinrath, JHEP 1111, 009 
(2011) jarXiv: 1108:06141 [hep-ph]]; S. F. Ge, D. A. Dicus and W. W. Repko, Phys. Rev. 
Lett. 108, 041801 (2012) [arXiv: 1108. 0964] [hep-ph]]: Riazuddin. [arXiv:1108.14M] [hep-ph]: 

F. Bazzocchi. larXiv: 1108.2497) [hep-ph]: T. Araki and C. -Q. Ceng, JHEP 1109, 139 (2011) 
[arXiv:1108!3T75] [hep-ph]]; S. Antusch, S. F. King, C. Luhn and M. Spinrath, Nucl. Phys. 
B 856, 328 (2012) jarXiv: 1108.42781 [hep-ph]]: A. Rashed and A. Datta, Phys. Rev. D 85, 
035019 (2012) [aFXiv:1109.2320l [hep-ph]]: P. O. Ludl, S. Morisi and E. Peinado, Nucl. Phys. 
B 857, 411 (2012) [arXiv:1109.3393l [hep-ph]]: G. Blankenburg and S. Morisi, JHEP 1201, 
016 (2012) jarXiv:1109.3396 [hep-ph]]: A. Aranda, C. Bonilla and A. D. Rojas, Phys. Rev. 
D 85, 036004 (2012) [a rXiv:1110.lT82] [hep-ph]]; G. -J. Ding, L. L. Everett and A. J. Stu- 
art, Nucl. Phys. B 857, 219 (2012) [arXiv:1110.1688] [hep-ph]]; D. Meloni, JHEP 1202, 
090 (2012) |arXiv:1110.5210l [hep-ph]]; S. Dev, S. Gupta, R. R. Gautam and L. Singh, 
Phys. Lett. B706, 168 (2011) jarXiv:1111.1300] [hep-ph]]; X. -G. He and S. K. Majee, 
larXiv: 1 111. 22931 [hep-ph] : A. Rashed. laFXiv: 1 1 1 1 .30721 [hep-ph] : R. d. A. Toorop, F. Feruglio 
and C. Hagedorn, Nucl. Phys. B 858, 437 (2012) [arXiv: 1112. 1 3401 [hep-ph] ] : S. F. King and 
C. Luhn. larXiv: 1 1 1 2 . 1959 1 [hep-ph] : S. Gupta, A. S. Joshipura and K. M. Patel, Phys. Rev. 
D 85, 031903 (2012) |arXiv:1112.6TT3l [hep-ph]]; A. Damanik, larXiv:1201. 27471 [hep-ph]; 

G. -J. Ding, larXiv:1201. 32791 [hep-ph] ; H. Ishimori and T. Kobavashi. la rXiv:1201. 34291 [hep- 
ph]; S. Dev, R. R. Gautam and L. Singh, Phys. Lett. B 708, 284 (2012) [a rXiv:1201. 37551 



19 



[hep-ph]]. 

[8] P. Minkowski, Phys. Lett. B67, 421 (1977); T. Yanagida in Workshop on Unified Theories, 
KEK Report 79-18, p. 95 (1979); M. Gell-Mann, P. Ramond and R. Slansky, Supergravity, 
p. 315; Amsterdam: North Holland (1979); S. L. Glashow, 1979 Cargese Summer Insti- 
tute on Quarks and Leptons, p. 687; New York: Plenum (1980); R. N. Mohapatra and 
G. Senjanovic, Phys. Rev. Lett. 44, 912 (1980). 

[9] K. S. Babu, R. N. Mohapatra, Phys. Rev. Lett. 70, 2845 (1993) |hep-ph/92092l5l . 
[10] G. Lazarides, Q. Shafi and C. Wetterich, Nucl. Phys. B 181, 287 (1981); J. Schechter and 
J. W. F. Valle, Phys. Rev. D 22, 2227 (1980); R. N. Mohapatra and G. Senjanovic, Phys. 
Rev. D 23, 165 (1981). 



[11] B. Bajc, G. Senjanovic and F. Vissani, |hep-ph/0110310| Phys. Rev. Lett. 90, 051802 (2003) 
[hep-ph/02102071 . 

[12] H. S. Goh, R. N. Mohapatra, S. P. Ng, Phys. Lett. B570, 215 (2003) |hep-ph/0303"055] ; 

H. S. Goh, R. N. Mohapatra, S. P. Ng, Phys. Rev. D68, 115008 (2003) |hep-ph/0308"l97] . 
[13] K.S. Babu, C. Macesanu, Phys. Rev. D72, 115003 (2005) |[hep-ph/0505"200l; S. Bertolini, 



M. Frigerio, M. Mahnsky, Phys. Rev. D70, 095002 (2004) [hep-ph/04061l7 |; S. Bertolini, 



T. Schwetz, M. Malinsky, Phys. Rev. D73, 115012 (2006) |; hep-ph/060 5006|; S. Bertolini, 
M. Mahnsky, Phys. Rev. D72, 055021 (2005) |hep-ph/0504241 1; G. Altarehi, G. Blanken- 
burg, JHEP 1103, 133 (2011) [arXiv: 1012.26971 [hep-ph]]; A. S. Joshipura, K. M. Patel, 
larXiv:1105.5"943 [hep-ph]. 
[14] C. S. Aulakh, B. Bajc, A. Melfo, G. Senjanovic, F. Vissani, Phys. Lett. B588, 196-202 
(2004) | hep-ph/0306242| ; T. Pukuyama, A. Ilakovac, T. Kikuchi, S. Meljanac, N. Okada, 
Phys. Rev. D72, 051701 (2005) |hep-ph/0412348] ; C. S. Aulakh and S. K. Garg, Nucl. 
Phys. B 857, 101 (2012) [ar Xiv:0807.0917l [heo-ph]]. 



[15] G. Altarelh and F. Feruglio, Rev. Mod. Phys. 82, 2701 (2010) [ arXiv:1002.0211 [hep-ph]]. 
[16] B. Dutta, Y. Mimura, R. N. Mohapatra, JHEP 1005, 034 (2010) jarXiv:0911.2242l [hep-ph]] . 
[17] S. F. King, C. Luhn, Nucl. Phys. B832, 414 (2010) |arXiv:0912.1344l [hep-ph]] : C. Hagedorn, 
S. F. King and C. Luhn, JHEP 1006, 048 (2010) [arXiv: 1003.42491 [hep-ph]]. 



20 



[18] See for instance, M.-C. Chen and K. T. Mahanthappa, Phys. Lett. B 652, 34 (2007) 
[arXiv:0705.0714l [hep-ph]]; C. Hagedorn, M. A. Schmidt and A. Y. Smirnov, Phys. Rev. 
D 79, 036002 (2009) [ arXiv:0811.2955] [hep-ph]]; H. Ishimori, Y. Shimizu and M. Tani- 
moto. Prog. Theor. Phys. 121, 769 (2009) |ar2OY:0812. 50311 [hep-ph]]; A. S. Joshipura, 
B. P. Kodrani, K. M. Patel, Phys. Rev. D79, 115017 (2009) |arXiv:0903.2161 [hep-ph]]; 
I. d. M. Varzielas, JHEP 1201, 097 (2012) |arXiv:1111.3952l [hep-ph]]. 

[19] P. F. Harrison, D. H. Perkins and W. G. Scott, Phys. Lett. B530, 167 (2002) 
[hep-ph/ 0202074|. 

[20] F. James and M. Roos, Comput. Phys. Commun. 10, 343 (1975); 
|http : //seal . web . cern . ch/seal/MathLibs/5_10/Minuit2/html/| 

[21] T. Blazek, S. Raby and S. Pokorski, Phys. Rev. D 52, 4151 (1995) fhep-ph/9504364l . 

[22] B. D. Wright, |hep-ph/9404217t M. S. Carena, S. Dimopoulos, C. E. M. Wagner and S. Raby, 
Phys. Rev. D 52, 4133 (1995) |hep-ph/ 9503488) ; D. M. Pierce, J. A. Bagger, K. T. Matchev 
and R. -j. Zhang, Nucl. Phys. B 491, 3 (1997) [hep-ph/9606"2lT] ; S. Antusch and M. Spin- 
rath, Phys. Rev. D 78, 075020 (2008) [ arXiv:0804.07T7l [heo-ph]]. 

[23] C. R. Das, M. K. Parida, Eur. Phys. J. C20, 121-137 (2001) [hep-ph/00 10004] . 

[24] J. Abdallah et al. (DELPHI Collaboration), Eur. Phys. J. C 55, 525 (2008); G. Rodrigo, 
A. Santamaria, M. S. Bilenky, Phys. Rev. Lett. 79, 193 (1997) [hep-ph/ 9703358]. 

[25] H. E. Haber and R. Hempfling, Phys. Rev. Lett. 66, 1815 (1991); J. Ellis, G. Ridolfi and 
F. Zwirner, Phys. Lett. B257, 83 (1991); R. Barbieri, M. Frigeni and F. Caravaglios, Phys. 
Lett. B258, 167 (1991); Y. Okada, M. Yamaguchi and T. Yanagida, Phys. Lett. B262, 54 
(1991); M. Drees and M. M. Nojiri, Phys. Rev. D 45, 2482 (1992); P. H. Chankowski, S. 
Pokorski and J. Rosiek, Phys. Lett. B274, 191 (1992). 

[26] F. E. Paige, S. D. Protopopescu, H. Baer and X. Tata, |hep-ph/0312045 [ 

[27] G. Aad et al. [ATLAS Collaboration], a rXiv: 1202. 14081 [hep-ex]: S. Chatrchyan et al. [CMS 
Collaboration], arXiv:1202. 14881 [hep-ex]. 

[28] S. Chatrchyan et al. [CMS Collaboration], Phys. Rev. Lett. 107, 221804 (2011) 
|a rXiv: 1109.2352l [hep-ex]]. 

[29] P. Nath and P. Fileviez Perez, Phys. Rept. 441, 191 (2007) |hep-ph/0601023] . 



21 



[30] K. S. Babu and S. M. Barr, Phys. Rev. D 48, 5354 (1993) | hep-ph/9306242] . 
[31] B. Dutta, Y. Mimura and R. N. Mohapatra, Phys. Rev. Lett. 94, 091804 (2005) 
|hep-ph/0412105|. 

[32] K. Abe et al, arX iv: 1109.32 62 [hep-ex]; Talk by M. Miura, 2011 Workshop on Baryon k 
Lepton Number Violation, Tennessee, US, September 22-24, 2011. 



22 



